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Evolution of the Ionizing Radiation Background and Star Formation in the 
Aftermath of Cosmological Reionization 
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O '■ ABSTRACT 
(N ; 
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^ . The temporal evolution of the ionizing UV background radiation field at high redshift 

. provides a probe of the evolution of the early star formation rate. By comparing the 

^ i observed levels of absorption in the highest redshift quasar spectra to the predictions of 

I a hydrodynamic simulation, we determine the evolution of the photoionization rate (F) 

I for neutral hydrogen in the intergalactic medium, over the redshift range 4.9 < z < 6.1. 

\jQ ' After accounting for sampling variance, we infer a sharp increase in F from z ~ 6.1 

■ to z ~ 5.8, probably implying reionization at this redshift. We find a decrease in F 

■ from z ~ 5.6 to 5.2, at 3a significance. This feature is a generic signature in the 
, aftermath of reionization, entirely consistent with the cosmological reionization process 

Q ! being completed at 2; ~ 6.1. 

fH I This generic feature is a result of a significant change in the star formation rate 

Qh| subsequent to the cosmological reionization. There is an abrupt rise of the temperature 

O ' of the intergalactic medium due to photo- heating, when it is reionized. This translates 

' to a correspondingly sudden jump in the Jeans mass and a dramatic suppression of 

d ■ gas accretion onto the most abundant (sub-galactic) halos at the epochs of interest. 

. The star formation rate suffers a temporary setback in the aftermath of reionization, 

^ , resulting in a temporary decrease in the amplitude of the ionizing radiation field. 



Subject headings: cosmology: theory — intergalactic medium — large-scale structure of 
universe — quasars: absorption lines 



1. Introduction 

The conventional wisdom based on the standard theory of structure formation says that the 
universe was reionized sometime in the redshift range z = 6 — 12 (Barkana &i Loeb 2001). The rela- 
tively large uncertain range in redshift refiects our imperfect knowledge of the density fiuctuations 
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on small scales, star formation processes (e.g., efficiency, IMF, etc.) and feedback processes at high 
redshift. The latest observations of high redshift quasars (e.g., Fan et al. 2001) are beginning to 
probe the lower bound of that window and suggestions have been made that we may be witnessing 
the end phase of the cosmological reionization process at z ~ 6, based solely on the appearance 
of a precipitous drop of transmitted flux at the rest-frame Lya wavelength near z ~ 6 in a single 
quasar spectrum (Becker et al. 2001, hereafter BOl; Barkana 2001). 

In this Letter we present an analysis of the ionizing background radiation field in the redshift 
range z = 4.9 — 6.1, using the absorption measurements in BOl. Combining with previous data at 
lower redshift, we find that, coming from high redshift, the ionizing radiation intensity displays a 
sharp rise at 2; ~ 6 peaking at 2; = 5.6, a significant downturn from z ~ 5.6 to 5.2 by a factor of 
~ 0.6, and subsequently a consistent ascent from z = 5.0 to z = 2.4 by a factor of ~ 4. All these 
three features are consistent with a picture that the cosmological reionization was near completion 
at 2; ~ 6. 

The first feature (i.e., the initial sharp rise) has been predicted by several authors previously 
(Cen & Ostriker 1993; Gnedin 2000a; Miralda^Escude, Haehnelt, & Rees 2001). At present, the 
primary uncertainty for its determination observationally is the possibility that the high level of 
absorption in the single observed quasar at z = 6.28 is some kind of anomaly (we show that it 
cannot be a simple statistical fluctuation in the absorption level). Additional observed quasars at 
z > 6.3 will be critical in this regard. Here we focus our attention on the relatively more robust 
measurements of the ionizing radiation field at 2; < 6 and show that the observed pause in the rise 
of the amplitude of the ionizing radiation field from z ~ 5.6 to 5.0 could be due to suppression of 
star formation following reionization at 2 ~ 6. 

2. Reionization and the Evolution of the Ionizing Background 

Assuming the model in which the IGM is almost completely ionized and the radiation back- 
ground is uniform, we infer the evolution of the H i ionization rate by requiring that the mean flux 

decrement in simulated spectra matches the observed values in the upper panel of Figure 2 of BOl. 
This method was pioneered by Ranch et al. (1997) and extended in McDonald et al. (2000,2001) 
and McDonald &: Miralda-Escude (2001). We apply the standard procedure using a hydrodynamic 
simulation of a flat universe dominated by a cosmological constant and cold dark matter (ACDM) , 
with CDM density $7^ = 0.3, Hubble parameter h = 0.67 {Hq = 100 /tkms^^/Mpc), baryon den- 
sity fi;, = 0.035, power spectrum normalization erg = 0.9, and large-scale primordial power spectrum 
slope n = 1. The simulation is Eulerian, with box size 25 h'^ Mpc divided into 768^ cells (see Cen 
et al. 2001 for a more complete description). We use outputs from the simulation at z = 5 and 
z = 6, interpolating linearly between them when necessary. We start by assuming the model in 
the simulation correctly represents the universe, and then consider whether or not our results are 
self-consistent. 
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Figure 1 shows the results for the evolution of r_i2(-2;), the hydrogen ionization rate in units 
of 10^^^ s^^. The open squares are the old results from McDonald & Miralda-Escude (2001). The 
filled squares arc the new results, leaving out the highest redshift point, which we analyze differently 
than the others because the detected transmission is not significant and the detector noise needs to 
be taken into account carefully (this last point is represented by the triangle with 1 and 2 a upper 
limit bars). The equally spaced bins in redshift are chosen to coincide with the majority of the 
points in Figure 2 of BOl. As discussed below, the results have been rescaled to apply for the value 
of the baryon density favored by big bang nucleosynthesis, rifo/i^ = 0.02 (Buries &; Tytler 1998), 
instead of the value in the simulation, Q,}jh? = 0.0157. 

In detail, we determine a value of r_i2 and its error for each BOl point using the simulation as 
follows: First, we construct the optical depth field, t{v) (i.e., the spectrum is F{v) = exp[— r(t;)]) 
along a large number of lines of sight through the simulation cube (256^ for each face), normalized 
using an arbitrary value of r_i2 (i.e., r oc for low density gas in ionization equilibrium). We 
then randomly choose many sets of five 12.5 hT^ Mpc-long segments of spectrum, so that each set 
forms a combined spectrum of the length of our redshift bins, Az = 0.167 (note that this length 
depends slightly on cosmology and redshift, i.e., the spectra should be shorter by 26% at z = 6 in 
our model, corresponding to a 13% increase in the estimated error bars, which we account for in 
the plot). For each 62.5 Mpc chunk, we determine the value of r_i2 that produces the observed 
value of the mean transmitted flux fraction, F{z) \T{z) in the notation of BOl]. The mean of the 
r_i2 values for all the chunks is the inferred value of r_i2 for this observed data point, and the 
dispersion is the estimated error. Most redshift bins in BOl have several measured values of F in 
them, so we combine the inferred values of r_i2 within the bin using an error weighted average 
which takes into account fractional overlap when necessary (because of the overlapping points, our 
error bars are correlated, in the sense that point-to-point differences are more significant than the 
error bars seem to indicate). 

To test that the treatment of 12.5 Mpc segments as independent does not underestimate 
the error by ignoring correlation along the lines of sight on larger scales, we repeat this procedure 
using 6.25 Mpc segments, and find almost identical error bars. 

We now turn to the highest redshift point in BOl. The given value and error are F = 0.0038 ± 
0.0026; however, BOl state that the sky subtraction has an uncertainty of order la, so we use 
a = 0.0052 when interpreting this point. BOl also give the measured transmitted flux fraction in 

the region of the spectrum corresponding to Ly/3 absorption by gas at this redshift, F = — 0.00024ib 
0.0024, and point out that the implied constraint on the Lya optical depth is actually stronger than 
the direct one. We assume the same sky subtraction error for this point as for the Lya region. 

We start by analyzing the Lya region. We use the method described above to create chunks 
of spectra of the appropriate length. For each value of r_i2 we compute F for many different 
chunks, and compute the mean likelihood of producing F = 0.0038 given the error a = 0.0052. The 
relative likelihood function peaks at r_i2 = 0.059, but does not decrease significantly for arbitrarily 
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Fig. 1. — Inferred ionization rate as a function of redshift. The open squares are from McDonald & 
Miralda-Escude (2001) while the filled squares and triangle are the new results in this paper. The 
triangle at the highest redshift is a non-detection with the error bars indicating 1 and 2a upper 
limits. Changing the cosmological parameters can cause coherent shifts in the overall amplitude 
of the points (not accounted for in the error bars), as can changing the temperature of the gas; 
however, none of these should cause variation in a small redshift interval, except if the temperature 
changed suddenly. 
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small r_i2. The likelihood has dropped by factors of 0.61 and 0.14 at r_i2 = 0.096 and 0.135, 
respectively, so we call these the 1 and 2(7 upper limits. 

We analyze the Ly/3 region in two ways. First we repeat the above procedure assuming a 
constant factor 0.13 suppression of the continuum by Lya absorption corresponding to gas at 
z ~ 5 (sec the lowest redshift point in BOl). The inferred 1 and la upper limits on r_i2 are 0.032 
and 0.056. These constraints are stronger than the ones from the Lycv region; however, they arc not 
as realistic as possible because they ignore fluctuations in the Lycc absorption in the Ly/? region, 
which might either weaken the constraints, by completely obscuring substantial portions of the 
Ly/3 region, or strengthen the constraints, by providing regions of significantly less than average 
absorption where residual transmission through the Ly/? forest could be observed. We simulate 
this effect by creating Lya absorption spectra to serve as the continua for the Ly/3 region, using 
the 2 = 5 simulation output and our inferred value of r_i2 from the lowest redshift BOl point. 
This procedure results in essentially identical constraints to the assumption of a flat continuum: 
r_i2 < 0.032 and 0.062 at 1 and 2cr, respectively. 

It may be possible to obtain even stronger constraints by using Ly7, as suggested by BOl 
(X. Fan, private communication), or by considering the limits imposed by individual pixel values 
(D. Weinberg, private communication). However, if reionization is on-going or just completed at 
z = 6.1, the detailed constraints on r_i2 obtained from any kind of analysis based on a simulation 
in which reionization happened much earlier probably are not very meaningful (r_i2 is not well 
defined if the radiation background is inhomogeneous) . Our formal constraints on r_i2 near z = 6 
should be taken simply as an indication that something dramatic is happening at this redshift and 
we think it likely indicates the end of the cosmological reionization process. 

The overall amplitude of our F results (but not the evolution over a short redshift interval) is 
sensitive to the cosmological model and gas temperature-density relation in the simulation, as one 
can see from the equation for optical depth at a point in real space (i.e., ignoring peculiar velocities 
and thermal broadening) with baryon overdensity A = pb/pb and gas temperature T: 

H{z) F ' ^ > 

where H{z) is the Hubble parameter at redshift z [H{z) ~ 100 h Q^rn^ {1+z)'^/'^ at high z in a ACDM 
universe] . Sensitivity of the inferred value of F to the power spectrum of density perturbations enters 
through the term, as well as the peculiar velocities, and additional sensitivity to temperature 
can arise from the thermal broadening. The effect of changing ri^^^ or H{z) on the inferred value 
of F is easy to correct for, but the effect of changing the temperature-density relation or the input 
power spectrum generally must be investigated by running additional simulations. We do this using 
approximate Hydro-PM simulations (hereafter HPM, Gnedin & Hui 1998), after checking that they 
give similar results to our fully hydrodynamic simulation (to ~ 10% in F) when performed using the 
same model and initial conditions. First, we analyze HPM simulations with a variety of different 
temperature-density relations and find that the thermal broadening is not very important to the 
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derivation of F, so the dependence F oc T~^-' is accurate (to ~ ±0.1 in the power law index), as 
long as the temperature is specified at the most relevant density for the transmission features that 
determine F, which we find to be A 0.35. The temperature in our hydrodynamic simulation 
at this density is 5400 K, but in the aftermath of reionization the temperature of the gas would 
be higher than this (e.g., Miralda-Escude & Rees 1994), requiring a correction. The amplitude of 
the initial density fluctuations in our hydrodynamic simulation is higher by a factor ~ 1.5 on the 
scale of the Lyct forest than the measurement of McDonald ct al. (2000), which was independently 
confirmed by Croft et al. (2001). By varying the amplitude. A, in HPM simulations, we determine 
that F oc A~^'^ (roughly). The required ionizing background decreases when A increases because 
the voids where transmission is observed become deeper (i.e., A becomes smaller). In the end, 
the corrections for power spectrum amplitude and temperature may roughly cancel (if we use 
T ~ 10000 K), so the points in Figure 1 should be reasonably accurate; however, we emphasize 
that the temperature at this redshift is unknown, so significant uncertainty in F remains. 

The reader may be wondering at this point "could the inferred decrease in F_i2 from z ~ 5.6 
to z ~ 5.2 possibly be real"? It appears to be a 3a effect and we think it could be real. The same 
trend is clearly visible to the eye in the bottom panel of Figure 2 of BOl, once we consider that 
the value of their Tg// is generally expected to increase with time. The signal is almost equally 
strong in the spectra of two different quasars. These two are at almost identical redshift so a 
continuum-related effect might be the same in both; however, the observed ^ 30% increase in 
transmission fraction between 2; ~ 5.2 and 2; ~ 5.6 would require an increase in the quasar 
continuum. Finally, we have every reason to believe that our statistical error calculation is correct. 
The scatter in the inferred values of F_i2 for different BOl points in each redshift bin is in perfect 
agreement with our estimated error bars (e.g., = 13-0 for 10 degrees of freedom when we treat 
the three average values of F_i2 for the bins in the range 5.2 < z < 5.6 as fitting parameters to 
be estimated from the individual points with error bars). Note that, while it may look small to 
the eye, the fractional error on the z = 5.25 point is actually larger than the fractional errors on 
the following two points. Purely adiabatic temperature evolution oc (1 + z)"^, which constitutes the 
steepest possible evolution in temperature (by ignoring the subsequent photo- heating), instead of 
the essentially constant temperature in the simulation, would only cause a 7.5% relative change in 
the values of the Fs. 

We think that the sharp rise of the radiation intensity at z ~ 6.1 may be a fairly strong 
case for the completion of cosmological reionization at that epoch. However, note that a value of 
F_i2 ~ 0.03 is not inconsistent in principle with a fully ionized IGM, since the neutral fraction 
near the mean density is still only ~ 0.0007. One is hampered by the high optical depth of an 
IGM that is still mostly ionized. It will be extremely valuable to look for additional signatures of 
cosmological reionization in its aftermath, where more flux is transmitted. We will therefore focus 
our discussion on the radiation field at z < 6. 
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3. Discussion 

In the standard hierarchical structure formation theory smaller structures began to form earlier 
than larger structures, simply because the amplitude of density fluctuations is a decreasing function 
of scale and in linear theory density fluctuations on all scales grow at the same rate. Under this 
general picture the universe is conventionally thought to be reionized mostly by photons from stellar 
systems more massive than 10^ M0, where atomic line cooling provides an efficient energy sink for 
gas that is collected and heated during gravitational collapse of halos (Haiman, Rees &: Loeb 1997; 
Gnedin & Ostriker 1997). Less massive systems, after having produced an insufficient amount of 
stars (Pop III) to reionize the universe at an earlier epoch, can no longer form stars due to lack of 
cooling processes (Haiman, Thoul, & Loeb 1996; Haiman, Rees & Loeb 1997; Tegmark et al. 1997). 

After reionization, the temperature of the intergalactic gas is raised to > lO^K (Cen & Ostriker 
1993; Miralda-Escude & Rees 1994; Gnedin & Ostriker 1997). The subsequent effect of suppression 
of gas accretion onto small halos has been noted by many authors (Efstathiou 1992; Thoul & 
Weinberg 1996; Quinn, Katz, & Efstathiou 1996; Kepner, Babul, & Spergel 1997; Navarro &: 
Steinmetz 1997; Kitayama &: Ikeuchi 2000; Gnedin 2000b). The extent of this effect somewhat 
varies among the studies and it is clear that our current treatment of the reionization process is 
far from ideal. Obviously, a sufficiently adequate study would require simulations that include an 
accurate treatment of three-dimensional radiative transfer process and have high enough spatial 
and mass resolution. We do not have simulations of this caliber at present and attempt to only 
make a semi-quantitative calculation to illustrate a consequence of reionization on star formation 
and the background radiation field, which we believe captures a primary feature of this process in 
a simple and intuitive fashion. 

Figure 2 shows the nonlinear mass and linear Jeans mass Mj 

cans ^ 

function of rcdshift. 

The nonlinear mass is the mass within a top-hat window at which the linear rms density 
fluctuation is unity, and Mjeans = 1-5 x 10^^{j^f/'^{l + z)^3/2q^i/2^-i jyj^ (Peebles 1993), 
where T is the gas temperature and CIm is the present matter density parameter. The assumption 
of instantaneous reionization is just for the sake of plotting convenience and is not required to make 
our point. In a real universe the reionization phase is complicated and is thought to progress on 
a time scale of a Hubble time, during which the mean (volume weighted) radiation field builds up 
slowly up to a value of approximately 10~^^ erg/cm^ /hz/sec/sr at Lyman limit, followed by a brief 
phase, when the majority of the baryons are ionized and a sudden jump in the amplitude of the 
mean radiation field intensity at Lyman limit to 10~^^ — 10"^^ erg/cm^ /hz/sec/sr occurs within a 
redshift interval of a fraction of unity (Miralda-Escude et al. 2000; Gnedin 2000a). We simply call 
this epoch of a sudden rise in radiation field the reionization epoch, marking the completion of the 
reionization process. 

We see, in Figure 2, that gas is able to accrete onto halos with mass greater than the Jeans 
mass ~ 10^ M0 prior to z ~ 6. Stars formed inside halos with M = 10^ — 5 x 10^ Mq are 
presumably primarily responsible for producing most of the ionizing photons immediately prior to 
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Fig. 2. — Nonlinear mass and linear Jeans mass Mjeans as a function of redshift, with the 
simplified assumption that the universe was instantaneously ionized at z = 6. The intergalactic 
medium is assumed to have the same temperature as the cosmic microwave background at z > 6 
and photo-heated to 1.5 x lO^Kelvin at z < 6. A ACDM model with n = 0.3, A = 0.7, n = 1 and 
CTg = 0.8 is assumed in order to compute M^i. 
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z ~ 6, where atomic line cooling provides an efficient cooling process. Right after reionization, 

there is a dramatic inversion of the situation from z ~ 6 to z ~ 4.6, where Jeans mass surpasses 
the nonUnear mass. Such a switch has some interesting consequences: gas accretion onto the most 
abundant halos (about the nonlinear mass) is suddenly suppressed, simply because the ambient gas 
has a temperature that is higher than the virial temperature of the typical halos and is incapable 
of cooling to a lower temperature. Gas accretion and thus star formation may be brought to a halt 
in typical halos (larger halos with higher virial temperature may still be able to form stars with 
a somewhat reduced rate due to the same effect). Below z ^ 4.6, nonlinear mass again becomes 
larger than the Jeans mass and gas is once again capable of accreting onto typical halos with mass 
> 4 X 10^ Mq, and star formation rate bounces back. We note that the cosmological model adopted 
may not be correct in fine details and is used only to illustrate the process. Therefore, the exact 
crossing redshift of the two curves at z = 4.6 should not be taken at its face value, but rather as a 
generic feature. 

Let us now relate the star formation rate to the ionizing background radiation field. Our 
treatment is intended to be simple but we think it contains the essence of the relevant physical 
processes. For gas immersed in a meta-galactic ionizing radiation background, locally hydrogen 
photoionization is balanced by hydrogen recombination (to levels n = 2 or higher): 

Vnm = nHnnea{T), (2) 

where T is the photoionization rate (see Figure 1), nni, uhii and rig are the neutral hydrogen, ion- 
ized hydrogen, and electron number densities, respectively, and a{T) is the hydrogen recombination 
rate. Averaging spatially both sides of equation (1) gives 

C<nH>^ a{T) 
< riHi >= , (3) 

where C is the clumping factor of all gas participating in receiving ionizing photons (i.e., excluding 
gas in optically shielded regions) . For simplicity we have assumed that the gas is composed entirely 

of hydrogen. Since the neutral hydrogen fraction is much less than unity, we have simply replaced 
both the ionized hydrogen and electron number densities by the total hydrogen number density 
hh- The mean free path of an ionizing photon in comoving length units is then 

A = ^±^^ ,,/" + -l . (4) 

where an is the hydrogen photoionization cross-section. Using a = 4 x 10~-^^cm^sec~''", T = 
1.5 X 10~^^ sec~^ (see Figure 1) and an = 2 x 10~^* cm^ we obtain A in comoving megaparsecs: 

X = 74C{z)-' (^^^ (^) ^°™i"g^P^- (5) 

Evidently, at z ~ 6, the mean free path of an ionizing photon is much smaller than the Hubble 
radius and radiative processes may be treated "locally" , by ignoring cosmological effects. But the 
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mean free path is much greater than the typical separation between ionizing sources 1 comoving 
Mpc). In this case wc may relate the mean comoving specific emissivity of ionizing radiation Cjy 
(in units of erg/sec/hz/comoving cm^) to the mean ionizing radiation intensity Ji^ (in units of 
erg/sec/hz/sr/cm^) by 

= 8.9 X l026c(z)-ie,(^)-2(i±^)-3. (6) 

The evolution of C{z) is complicated and currently undetermined without detailed simulations 
that ionized the universe at 2; ~ 6. However, simple estimates may be made to gauge the range 
of variation. In the linear regime we have C{z) cx (1 + z)""^ and in the extreme nonlinear regime 
C{z) oc (1 + z)^; the actual situation is likely to be between these two limiting cases, parameterized 
as C{z) oc (1 + z)~'^~^^ with j3 > 0. Inserting C{z) into equation (5) yields 

oc e^(l + z)-^-^ (7) 

with /3 > 0. Clearly, in order for Ji, to increase with redshift, as indicated by observations in the 
redshift range z = 5.2 — 5.6 as shown in Figure 1, the comoving specific emissivity or comoving 
star formation rate has to increase with redshift. Taking the values shown in Figure 1 in the 
range z = 5.6 — 5.2, we obtain oc (1 + z)^-^^^-^. (This fit used the 3rd through 5th filled 
squares in Figure 1. We obtain a flatter slope, 5.3 ± 2.1, with an acceptable if we add the 2nd 
point, at z ~ 5.1.) This requires a decrease in star formation rate from z = 5.6 to z = 5.2 as 
e oc (1 + z')^0-^^^-^-l^ _ Our analysis is in qualitative agreement with Barkana Sz Loeb (2000), who 
pointed out that the decrease in star formation subsequent to cosmological reionization may be 
detectable in the evolution of number counts of faint galaxies by NGST. It may be that we have 
detected this signature in the evolution Ji, or F. 



4. Conclusions 

In hierarchical structure formation theory the evolution of the meta-galactic ionizing radiation 
background may be characterized by five distinct phases: 1) from z ~ 15 — 20 to Zri (which marks 
the end of reionization) the radiation field builds up slowly before individual ionizing HII regions 
overlap ( "pre-overlap period"; Gnedin 2000a) 2) in a brief period up to Zri the majority of baryons 
are ionized in a short time scale and the radiation field jumps up by about two to three orders 
of magnitude, completing the reionization process ("reionization epoch"); 3) from Zri to zjeans 
the radiation field pulls back significantly due to a decrease in star formation rate as a result of 
reionization (we call this "aftermath of reionization" ) ; 4) zj^ans to z^eafc ~ 1 — 2 the radiation field 
rises steadily with contributions from both stars in larger galaxies and quasars; 5) Zpgafe 1 — 2 
to z = the radiation field drops off sharply due to the combined effect of decreasing star /quasar 
formation rate and cosmological effects. Consistent with theoretical calculations, observations have 
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confirmed phases (4) (e.g., McDonald & Miralda-Escude 2001; shown in Figure 1 of this paper) 
and (5) (e.g., Shull et al. 1999). 

In this paper we have interpreted the absorption in the spectra of the latest high redshift 
quasars in the redshift range z = 4.9 — 6.1, searching for salient features pertaining to the "reion- 
ization epoch" (phase 2) and the "aftermath of reionization" (phase 3). Quite intriguingly, we 
appear to be seeing the "aftermath of reionization" and the end of the "reionization epoch" in 
these spectra, both of which seem to point to a cosmological reionization redshift of Zri ~ 6. Specif- 
ically, a predicted decrease of star formation rate in the aftermath of the reionization is indicated 
by the observations. While a few more high redshift quasars are required to check that the com- 
plete absorption at z = 6.1 is not some kind of anomaly (e.g., a rare leftover neutral patch in an 
almost completely ionized universe), it appears that our conclusion with regard to the "aftermath 
of reionization" is more secure due to its independence of any single spectrum. It is still very urgent 
to enlarge the observational data sets to greatly firm up the conclusion and sharpen up the rate 
of evolution. If confirmed by future observations, this may be indication that our current under- 
standing of galaxy formation may be approximately valid up to redshift z ~ 6. Furthermore, as 
suggested by Barkana & Loeb (2000), it will be important to independently detect this "aftermath 
of reionization" signature in the evolution of number counts of faint galaxies by NGST. 

Since we appear to be observing the period directly after reionization, it becomes more impor- 
tant to perform dedicated simulations that achieve reionization at the appropriate redshift z ~ 6. 
This paper should be understood as an attempt to guide that effort by making some preliminary 
comparisons between data and an existing simulation in which reionization occurred somewhat ear- 
lier. Such future simulations are quite demanding, because not only high resolutions are required 
but also an accurate treatment of three-dimensional radiative transfer must be implemented. 

This research is supported in part by grants AST93-18185 and ASC97-40300. We thank 
Xiaohui Fan, Vijay Narayanan and David Weinberg for useful discussion, and Nick Gnedin for his 
HPM code. 
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